Suppression of Dimer Correlations in the Two-Dimensional J1-J2 Heisenberg Model: 

an Exact Diagonalization Study 
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(Dated: February 2, 2008) 

We present an exact diagonalization study of the ground state of the spin- half J\ — Ji model. 
Dimer correlation functions and the susceptibility associated to the breaking of the translational 
invariance are calculated for the 4x4 and the 6x6 clusters. These results - especially when 
compared to the one dimensional case, where the occurrence of a dimerized phase for large enough 
frustration is well established - suggest either a homogeneous spin liquid or, possibly, a dimerized 
state with a rather small order parameter. 

PACS numbers: 75.10.Jm, 71.27.+a, 74.20.Mn 



In quantum systems, frustration, i.e., the impossibil- 
ity to satisfy simultaneously every pairwise interaction, 
induces unconventional properties down to zero tempera- 
ture (T = 0). In particular it plays a very important role 
in quantum antiferromagnets where the combined effect 
of competing interactions and zero-point motion can give 
rise to non-magnetic ground states which lack the clas- 
sical Neel long-range order. To this respect the most 
intriguing situation appears in two spatial dimensions 
where at T — both the magnetic and non-magnetic 
ground states are allowed by the Mermin- Wagner theo- 
rem, and quantum phase transitions from an ordered to a 
disordered magnetic phase can be triggered by increasing 
the frustration^ 

One of the simplest frustrated spin systems in two di- 
mensions is the so-called J1 — J2 Heisenberg models Here 
the antifcrromagnctic coupling (Jj > 0) between nearest 
neighbor (n.n.) spins is frustrated by the presence of 
a next-nearest neighbor (n.n.n.) superexchange interac- 
tion (J 2 > 0): 



n = Jr s< • s,- + J 2 



Si • S i 



(i) 



where Sj are spin-half operators on a periodic N— site 
square lattice. Recently, the interest in this model has 
been boosted by the synthesis of three vanadates com- 
pounds (Li 2 VOSi0 4 , Li 2 VOGe0 4 , and VOM0O4) whose 
relevant magnetic interactions involve nearest and next- 
nearest spin- 1/2 F 4+ ions on weakly coupled stacked 
planes £ 

Despite the simplicity of the Hamiltonian the 
zero-temperature phase diagram of the spin- half J1 — J2 
model has been much debated in the last 15 years. For 
J2/J1 *C 0.5, an antifcrromagnctic (AF) Neel order is 
expected, while in the opposite limit, J2/J1 3> 0.5 the 
ground state should be in an AF collinear phase, with 
the spins ferromagnetically aligned in one direction and 
antifcrromagnetically in the other, corresponding to mag- 
netic wave vectors Q = (n, 0) or Q = (0,tt). The most 
interesting part of the phase diagram is the one around 



the fully frustrated point J2/J1 = 0.5. In fact, although 
the existence of a gapped (non-magnetic) phase is at 
present rather likely for 0.4 < J2/J1 < 0.55^^ its char- 
acterization is one of the most intriguing puzzles of the 
physics of strongly correlated systems. In particular, a 
debated issue is whether the ground state is a homo- 
geneous spin liquid, as originally proposed by Anderson 
and FazckasAi The other possibility, following Read and 
Sachdev's large- iV expansion^ is a spontaneously dimer- 
ized phase which breaks the translation invariance of the 
Hamiltonian or - more in general - one of the lattice 
symmetries i^iiMi 

The main reason for the lack of agreement on this 
delicate issue is due to the fact that most of the ex- 
isting results have been obtained with approximated 
numerical and analytical techniques based on reference 
states explicitly breaking some symmetry of the Hamil- 
tonianj 4 ! 5 ! 9 ! 10 : 11 so that it is very difficult to put on a 
solid ground any conclusion on the actual ground-state 
correlations. In this context, the unbiased results pro- 
vided by the exact diagonalization (ED) of small clusters 
can give a valuable insight into the T = properties of 
the models In particular, on the 4x4 cluster, Poilblanc 
and collaborators^ showed that the dimer correlations 
are enhanced in the parameter range 0.5 < J2/J1 < 0.6, 
concluding that the dimer state is a good candidate for 
describing the physics of the J1 — J2 model in the highly- 
frustrated region. This claim is questionable because of 
the smallness of the 4x4 lattice, where it is very diffi- 
cult to extract the important long-range behavior of the 
correlation functions. In this paper, we enlarge the clus- 
ter up to N = 36 sites and present a more systematic 
analysis based on i) the calculation of the dimer suscep- 
tibility in the presence of an explicit symmetry-breaking 
term and ii) the calculation of a static dimer correla- 
tion function. In both cases, the direct comparison with 
the one-dimensional chain, where the dimerization is well 
established for large enough frustration^ casts doubts 
on the existence of a dimerized two-dimensional ground 
state in the non-magnetic regime. 
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FIG. 1: ED calculation of the susceptibility associated to 
the operator Ot of Eq. P). Left panels: results for the one- 
dimensional chain. Bottom: Xt(S)/N vs S 2 for (from below) 
4, 6, and 8 sites. J2/J1 =0.1 (empty symbols), 0.4 (full sym- 
bols). Top: size scaling of xt/N. Right panels: same quan- 
tities for the two-dimensional case for J2/J1 = 0.5 (empty 
symbols) and 0.55 (full symbols). 



The occurrence of some kind of dimerized order can 
be detected by calculating the response of the system to 
operators breaking the lattice symmetries. This can be 
done by adding to the Hamiltonian a term —SO, where 
= Yli Oi is an (extensive) hermitian operator that 
breaks some symmetry of TL and 8 is a (small) parameter. 
In fact, if true long-range order in the dimer correlations 
exists in the thermodynamic ground state, the finite- 
size susceptibility \o = 2(V>o|0(i?o - H)^ 1 0\'ipo)/NJi 
has to diverge with the system size. In particular, it 
can be shown that it is bounded from below by the 
system volume squared, \o > const x p 4 N 2 , where 



V = Y ('0o|O 2 |V'o)/W 2 is the order parameter^S*^ A 
spontaneously broken translation symmetry can be de- 
tected by studying the response of the system to the per- 
turbation 50t, with 



j+X ) 



(2) 



where Rj indicates the coordinates of the site j and 
x = (1,0). This operator preserves the SU(2) symmetry 
of the J1 — J2 Hamiltonian but breaks the translational 
invariance with momentum Q = (n, 0) and the x «-> y 
symmetry (or equivalently the it/2 rotation symmetry). 

By using a numerical technique, like the Lanczos 
method, the susceptibility xt — —d 2 e(5)/dS 2 \s= can be 
calculated with only energy measurements by computing 
the ground-state energy per site (in unit of J\) in pres- 




FIG. 2: ED results of the dimer-dimer correlation functions 
for the 32-site one- dimensional chain with J2/J1 =0.1 (empty 
dots), 0.4 (full dots), and 0.5 (stars). Notice however that at 
the exactly solvable point J2/J1 = 0.5, |A Xia; (r)j = 1/64, 
independent on the distance for any non-overlapping singlets. 



ence of the perturbation, e(S), for few values of S and by 
estimating numerically the limit xt = lim,5_>o Xt(^) = 
— 2(e(6) — e )/S 2 . Notice that, in the presence of the 
perturbation, Eq. J5J), the Lanczos calculation is rather 
demanding. In fact, the symmetrized Hilbert space is in- 
creased by a factor of four with respect of the standard 
calculation on the J1 — J2 model ji£ and the resulting di- 
mension of the ground-state subspacc is 63.117.760 for 
N = 6 x 6. 

In Fig. n we present the ED results for the suscep- 
tibility associated to the operator Ot- For compari- 
son, the data for the one-dimensional J1—J2 model are 
also shown. In the latter case, it is well known that 
a transition form a quasi-ordered state, with power-law 
spin-spin correlation functions, to a dimerized phase oc- 
curs at J2/J1 — 0.241.M The response of the system to 
the perturbation is very different below and above the 
critical point [Fig. [TJa)-(b)]. This is particularly evi- 
dent by performing the size-scaling of Xt/N which di- 
verges for J2/J1 = 0.4, while saturates to a constant 
for J2/J1 =0.1 consistently with the known scaling di- 
mension of the spin-Peicrls perturbation in the Luttingcr 
regime, X = 1/2 [Fig. ^b)]»i£ In contrast in the two- 
dimensional case, for both J2/J1 = 0.5 and J2/J1 = 0.55, 
Xt/N decreases by increasing the lattice size [Fig.^c)- 
(d)] suggesting that a divergence of this quantity in the 
thermodynamic limit is unlikely, even if not ruled out 
by our small-size calculations. Indeed, Sushkov et alm& 
found a divergence of the same susceptibility xt, by us- 
ing a series expansion, that are free from size effects but 
are strongly biased by the approximate approach: the 
limited number of coefficients known in the infinite series 
and the difficulty of reconstructing a possible transla- 
tional invariant state starting from a dimerized reference 
one. Note that two-dimensional systems behave very 
differently from chains of comparable linear size, where 
Xt/N always increases with N in the dimerized phase 
[Fig. nja)-(b)]. This points toward the absence of dimer 
order parameter in the strong frustration regime. 

In order to investigate the possible occurrence of a 
spin-Peicrls phase we now analyze the dimer-dimer corre- 
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FIG. 3: ED calculation of the dimer-dimer correlations 
among pairs parallel to the y axis as a function of the Manhat- 
tan distance d for iV = 6 x 6. J2/J1 = 0.2 (empty dots), 0.45 
(empty squares), 0.55 (empty triangles), 0.6 (full triangles), 
0.65 (full squares), and 10 (full dots). The stars (crosses) in- 
dicate the results for the lowest d-wave (s-wave) excited state 
for J2/J1 = 0.65 (J2/J1 — 10). The large circles in the top 
panel are the values for J\ = 0. 

lation functions of the unperturbed (5 = 0) ground state: 

^n,v{r — t ) = {S*Sr +tl S*tS*t + J) — (S* S* +ji ) {S*i Sp +v ) , 

where \x and v are two unit vectors connecting the 
nearest-neighbor sites on the lattice. In presence of 
a broken spatial symmetry, the latter has to converge 
to a finite value for large distance, displaying also the 
typical staggered pattern shown in Fig. [3 for the one- 
dimensional models In this case, the different decay of 
the dimer correlations below and above the transition can 
be also recognized: oscillatory power-law in the Luttinger 
regime and constant amplitude oscillations in the dimer- 
ized phase. In Fig.[3]wc present a systematic study of the 
dimer-dimer correlations on the 6x6 cluster as a function 
of J2/J1. The lower panel is for J 2 /Ji = 0.2 (AF Necl 
phase), the middle one is for 0.5 < J2/J1 < 0.65 and the 
panel on the top is for J2/J1 = 10 (AF collinear phase). 

The dimer correlations do not seem to be much affected 
entering the non-magnetic region (J2/J1 > 0.4) and the 
large distance value of the correlations for J2/J1 — 0.5 is 
comparable to the one in the AF Neel phase, where the 
dimer correlations decay to zero with a power law. By in- 
creasing further the frustration, the dimer-dimer correla- 
tions increase reaching a maximum around J2/J1 — 0.65 
and then converging for large J2/J1 to a non-zero value, 
which is positive for parallel dimcrs and negative for or- 
thogonal ones (not shown). Notice, however, that for 
large enough frustration, A Vty (r) does not show a stag- 



gered behavior, which is inconsistent with the pattern 
predicted on the basis of a standard spin-Peierls state. 
On the other hand, the broken translation symmetry 
of the collinear phase, also expected in the large- J 2 re- 
gion, cannot be detected by the dimer-dimer correla- 
tion functions plotted in Fig. Instead, the fact that 
A M .j,(r) has different signs for parallel and orthogonal 
dimers is the signature of the breaking of the 7r/2 rota- 
tion symmetry, which also characterizes the AF collinear 
phase. In fact, for J 2 3> Ji, the system decouples in 
two unfrustrated Heisenberg models and A^ v (r — r') ~ 
{S'S^,){S' +ll Sp +v ) with the two pairs belonging to op- 
posite sublattices and (S*S*,) being the ground-state 
spin correlations of the Heisenberg model on a iV/2-site 
cluster (see the top panel of Fig- - Hence, in this limit, 
A Mi „(r - r>) ~ 8w\{$'fy)& +li fy+ v )\ with v , = 1 
and —1 for parallel and orthogonal dimers, respectively. 
Interestingly enough, for J 2 3> Ji, the large distance 
value of the dimer correlations will approach in the ther- 
modynamic limit A AIjly (r) ~ (m^) 4 /9, where m) ~ 0.307 
is the antiferromagnctic order parameter for the spin-half 
Heisenberg models 

The non-monotonic behavior of the dimer correlations, 
displaying a maximum around J2/J1 — 0.65 is likely to 
be a finite-size effect related to the finite-size gap between 
the two lowest singlet states. These states cross on the 6x 
6 for J2/J1 — 0.65, the ground-state symmetry changing 
from s-wave (below) to d-wave (above) H In fact, as we 
have checked in the AF collinear phase, the value of the 
dimer-dimer correlations on the s-wave singlet (crosses 
in the top panel of Fig. |3J) is larger than the value on the 
d-wave one, the two values converging to the same only 
when they become degenerate, i.e., at the level crossing 
on finite-sizes (stars in the middle panel of Fig. [2J, and 
in all the AF collinear phase in the thermodynamic limit. 

Following Refi^i, it is possible to obtain a finite-size 
estimate of the dimer order parameter from the long- 
distance behavior of the dimer-dimer correlations shown 
in Fig. □ 

D % d = 9 lim \(A y , v (r - y) - 2A y , y (r) + A y , v (r + y)\, 

\r\ — >oc 

where the factor 9 is required to take into account the 
three spin component of the order parameter. The be- 
havior of Dd as a function of J2 / J\ is reported in Fig.Qfa) 
for the 4x4 and 6x6 lattices. A remarkable decrease can 
be appreciated by increasing the lattice size, even in the 
strongly frustrated regime. Moreover, the maximum or- 
der parameter on the 6x6 lattice is considerably smaller 
than the corresponding one for the dimcrized chain, i.e., 
Dd = 0.75 for the Majumdar-Ghosh case (J2/J1 = 0.5). 
These facts should suggest that a sizable dimer order pa- 
rameter does not survive in the thermodynamic limit, 
leaving the possibility of a disordered spin liquid (or, at 
most, of a ground state with a very small dimer order). 

Finally, in analogy with the definition of Dd, it is pos- 
sible to define a dimer order parameter related to a ro- 
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FIG. 4: (a): ED calculation of the finite-size the dimer order 
parameter Dd as a function of J2/J1, for the 4x4 (full tri- 
angles) and the 6x6 (full circles) lattice, (b): The same for 

Do. 



tational symmetry breaking: 

D 2 a = 18 lim \(A y , y (r) - A v>x (r)\ 

\r\ — >oo 



Again, the behavior of D Q as a function of J2/J1 is re- 
ported in Fig. 0Jb) for the 4x4 and 6x6 lattices. A 
very sharp increase of D Q around J2/J1 ~ 0.6 confirms 
that a rotational symmetry breaking is plausible for large 
frustration, where the collinear phase is settled down^ 
For smaller values of J2/J1, instead, the ground state is 
rotational invariant, as previously found^ 

In conclusion, we have presented an exact diagonaliza- 
tion study of dimer susceptibility and dimer correlations 
of the spin-half J\ — J2 model. Although we are aware 
that important finite-size effects occur near phase tran- 
sitions, our analysis does not provide evidence in favor 
of a spin-Peierls ground state. By contrast, a homoge- 
neous spin-liquid ground state may be realized in the 
regime of strong frustration. This possibility has been 
recently supported by a variational calculation showing 
that a spin-liquid projected BCS wave function provides 
an almost exact variational ansatz of the ground-state 
in the non-magnetic phase of this frustrated system^ 
Work is in progress to corroborate these conclusions and 
to clarify the nature of the excitation spectrum of this 
unconventional state of matter. 
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